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For the determination of the 144 Sm-a optical potential we measured the angular distribution of 
144 Sm(a,a) 144 Sm scattering at the energy Ei^h = 20 MeV with high accuracy. Using the known 
systematics of a-nucleus optical potentials we are able to derive the 144 Sm-a optical potential at 
the astrophysically relevant energy -E c . m . = 9.5 MeV with very limited uncertainties. 
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I. INTRODUCTION 



In a detailed study of nucleosynthesis in Type II supernovae Woosley and Howard proposed the so-called 7-process, 
which is important for the production of the samarium isotopes 144 Sm and 146 Sm jl]. The existence of the 7-process 
was confirmed later Q]. Because of the a decay of 146 Sm (Ti/ 2 = 1-03 x 10 s y) today one can find correlations between 
the 142 Nd/ 144 Nd ratio and the Sm/Nd ratio in some meteorites || as a consequence of the 7-process. The nuclear 
reaction rates used in the network calculations of the 7-process are relatively uncertain. Especially, the reaction rate 
for the production of 144 Sm by the photodisintegration reaction 148 Gd(7,a) 144 Sm is uncertain by a factor of 10 [0. 

The reaction rate for the reaction 148 Gd(7,o) 144 Sm was derived from the 144 Sm(a,7) 148 Gd reaction cross section 
using a detailed-balance calculation (see e.g. ^H|). Usually, reaction rates are given at the relatively high temperatures 
of 2.5 < T 9 < 3.0, or the capture cross section at E c , m . — 9.5 MeV (corresponding to Tg = 2.8) is derived jllH-m]. 

Two ingredients enter into the Hauser-Feshbach calculation of the 144 Sm(a,7) 148 Gd capture cross section: transition 
probabilities and the nuclear level density. The transition probabilities were calculated using optical wave functions 
in an equivalent square well (ESW) potential (|^| , a Woods-Saxon (WS) potential j^J^] , and a folding potential HQ . 

Calculations with ESW potentials are very sensitive on a proper choice of the radius parameter R; two calculations 
using R = 8.01 fm (Ref. Jl, based on the ESW radius from Ref. j|) and R = 8.75 fm (Ref. §, based on Ref. §) differ 
by a factor of 10 for the 4 Sm(a,7) 148 Gd cross section. The WS potential and the folding potential in Rcfs. [?]] 
were taken from global parametrizations of a optical potentials; these results for the 144 Sm(a,7) 148 Gd cross section 
lie between the different ESW results. First experimental results on the 144 Sm(a,7) 148 Gd capture cross section at 
somewhat higher energies lie at the lower end of the different calculations fl0[| . 

The aim of this work is to determine the optical potential at the relevant energy E cm . = 9.5 MeV. In general, 
optical potentials can be derived from elastic scattering angular distributions. Recently, in a systematic study the 
energy and mass dependence of a-nucleus potentials was determined JO). In that work a scattering on 144 Sm was 
analyzed at higher energies. An extrapolation to astrophysically relevant energies is possible only with very limited 
accuracy. At the energy E c . m , — 9.5 MeV the 144 Sm(a,a) 144 Sm scattering cross section is given by almost pure 
Rutherford scattering because of the height of the Coulomb barrier of about 20 MeV. For a reliable determination 
of the optical potential one has to increase the energy in the scattering experiment. However, because of the energy 
dependence of the optical potential, the energy should be as close as possible to the astrophysically relevant energy 
Ec.m. = 9.5 MeV. As a compromise we measured the 144 Sm(a,a) 144 Sm angular distribution at £"i a b = 20 MeV. At 
this energy the influence of the nuclear potential on the angular distribution is measurable, even though it is small. 



* APART-fellow, present adress: Institut fur theoretische Physik, Universitat Basel, Klingelbergstr. 82, CH-4056 Basel, 
Switzerland 

^present address: Institut fur Strahlenphysik, Universitat Stuttgart, Allmandring 3, D-70569 Stuttgart, Germany 



1 



Especially the determination of the shape of the optical potential remains difficult even at this energy. For the real 
part this problem vanishes because its shape is given by a folding procedure, but the shape of the imaginary part has 
to be adjusted to the experimental angular distribution. For that reason the angular distribution has to be determined 
with very high accuracy in the full angular range. 

The increase of the real part of the optical potential at energies close to the Coulomb barrier is well-known |l2] . Close 
to the Coulomb barrier the number of open reaction channels changes strongly; this leads to a strong variation of the 
imaginary potential, and the strength of the real potential is coupled to the imaginary part by a dispersion relation. 
A further influence on the potential strength coming from antisymmetrization effects is indicated by microscopic 
calculations which were performed for light nuclei. This "threshold anomaly" (mainly so-called in heavy-ion scattering 
and fusion) was seen also in a scattering on many nuclei []il]-[l6f . 



II. EXPERIMENTAL SETUP AND DATA ANALYSIS 



The experiment was performed at the cyclotron laboratory at ATOMKI, Debrecen. We used the 78.8 cm diameter 
scattering chamber which is described in detail in Ref. |17|] . Here we will discuss only those properties which are 
important for our experiment. 



A. Targets and Scattering Chamber 

The samarium targets were produced by the reductive evaporation method (T^| at the target laboratory at ATOMKI 
directly before the beamtime to avoid the oxidation of the metallic samarium. A thin carbon foil (thickness d w 
20 ^g/cm 2 ) was used as backing. 

During the experiment we used one enriched 144 Sm and one natural Sm target, and one carbon backing without 
samarium layer. The enrichment in 144 Sm was (96.52±0.03)%. The thickness of the samarium targets was determined 
by the energy loss of the a particles in the samarium layer. We compared the energy of a particles elastically scattered 
on 12 C in the carbon backing at z?i a b = 162° using the pure carbon target as reference and both samarium targets. 
The stopping powers dE/dx of the a particles in samarium at the relevant energies E = 20 MeV (incident a) and E 
= 5.2 MeV (backward scattered a) were taken from Ref. The resulting thicknesses are d — 142 /ig/cm 2 ( 144 Sm) 
and d = 218 /ig/cm 2 (natura l Sm ) with uncertainties of about 10%. The pure carbon target was used also for the 
angular calibration (see Sect. II C| ). 



Additionally, two apertures were mounted on the target holder to check the beam position and the size of the 
beam spot directly at the position of the target. The smaller aperture had a width and height of 2 mm and 6 mm, 
respectively. This aperture was placed at the target position instead of the Sm target before and after each variation 
of the beam current. Because practically no current could be measured on this aperture the width of the beam spot 
was definitely smaller than 2 mm during the whole experiment, which is very important for the precise determination 
of the scattering angle. In contrast, the relatively poor det ermination of the height of the beam spot does not disturb 
the claimed precision of the scattering angle (see Sect. |ll C| ). Furthermore, the position of the beam on the target was 
continuously controlled by two monitor detectors. No e viden ce was found for a change of the position by determining 
the ratio of the count rates in both detectors (see Sect. |ll B ). 



B. Detectors 



For the measurement of the angular distribution we used 4 silicon surface-barrier detectors with an active area 
A = 50 mm 2 and thicknesses between d — 300 /jm and d — 1500 /im. The detectors were mounted on an upper 
and a lower turntable, which can be moved independently. On each turntable two detectors were mounted at an 
angular distance of 10°. Directly in front of the detectors apertures were placed with the dimensions 1.25 mm x 5.0 
mm (lower detectors) and 1.0 mm x 6.0 mm (upper detectors). Together with the distance from the center of the 
scattering chamber d = 195.6 mm (lower detectors) and d — 196.7 mm (upper detectors) this results in solid angles 
from Ail = 1.63 x 10 -4 to Ail = 1.55 x 10~ 4 . The ratios of the solid angles of the different detectors were determined 
by overlap measurements with an accuracy much better than 1%. 

Additionally, two detectors were mounted at the wall of the scattering chamber at a fixed angle of d = 15° (left and 
right side relative to the beam direction). The solid angles of these detectors are Ail — 8.10 x 10~ 6 . These detectors 
were used as monitor detectors during the whole experiment. 
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The signals from all detectors were processed using charge-sensitive preamplifiers (PA) , which were mounted directly 
at the scattering chamber. The output signal was further amplified by a main amplifier (MA). The bipolar output 
of the MA was used by a Timing Single Channel Analyzer (TSCA) to select signals with amplitudes between « 1 
V and 10 V, and the unipolar output of the MA was gated with the TSCA signal using a Linear Gate Stretcher 
(LGS). The LGS output was fed into an CAMAC ADC, and the ADC data were stored in a corresponding CAMAC 
Histogramming Memory module. The data acquisition was controlled by a standard PC with a CAMAC interface 
using the program MCMAIN |2C| . The dead time of the system was determined by test pulses, which were fed into 
the test input of each PA. It turned out that the deadtime was negligible (< 0.2%) except the runs at very forward 
angles. 

The achieved energy resolution was better than 0.5% corresponding AE < 100 keV at E a » 20 MeV. 



C. Angular Calibration 



Because of the strong angular dependence of the scattering cross section especially at forward angles the angular 
calibration was done very carefully by two kinematic methods. The carbon backing contained some hydrogen contam- 
ination. Therefore, we used the steep kinematics of 1 H(a,a) 1 H scattering at forward angles (10° < $i a b < 15°). We 
measured the energies of the a particles scattered from 1 H (note: two peaks with different a energies corre spon ding 
to two center-of-mass scattering angles can be found at one laboratory scattering angle, labelled ± in Eq. |2.l| ) and 
from 12 C (ground state and 2 + state at 4.44 MeV), and we determined the ratio 

± m _ -g Q ( 12 Cg. s .)-^( 1 H) 

from the experimentally measured energies and from a calculation of the reaction kinematics. The angular offset is 
given by the mean value of the differences in i? derived from all determined values of q^ xp and qf aic - The following 
results can be obtained from this procedure: A-i? ff sc t = —0.38° ±0.02° (lower detectors) and Az9 ff sc t = +0.32° ±0.02° 
(upper detectors). (Note that these offsets cannot be a consequence of a beam spot which is not exactly centered on 
the target; for the measured offset of about 0.35° the beam spot would have been outside the aperture with a width 



of 2 mm, which can be placed at the position of the target, see Sect. [I A.) The uncertainty of the adjustment of the 
angle is given by the standard deviation of each single measurement: A$ ad j USt < ±0.1° (all detectors). 

In a second step we measured a kinematic coincidence between elastically scattered a particles and the corresponding 
12 C recoil nuclei. One detector was placed at #i a b,a = 77° (lower detector, left side relative to the beam axis), and 
the signals from elastically scattered a particles on 12 C were selected by an additional TSCA. This TSCA output 
was used as gate for the signals from another detector which was moved around the corresponding 12 C recoil angle 
$iab, recoil = 42° (upper detector, right side). The maximum recoil count rate was found almost exactly at the expected 
angle (see Fig. [l]) . 



D. Experimental Procedure, Uncertainties 

With this setup we measured spectra from the 144 Sm and the natural Sm target at angles from 15° to 172° in 
steps of 1° (■& < 140°) and 2° (i? > 140°). Two typical spectra measured at forward ($ = 25°) and at backward 
(i) = 130°) angles are shown in Fig. 0. The measuring times were between some seconds (forward angles) and several 
hours (backward angles), and the corresponding beam currents were between 30 nA and 600 nA 4 He 2+ ions. The 
beam was stopped in a Faraday cup roughly 2 m behind the scattering chamber, and the current was measured by a 
current integrator. 

For each run the scattering cross section was determined relative to the monitor count rate: 



is\ — / — \ is — . : : v " 1 . — Mon (o 2) 

\dnJ y ' WfJ/Mon ' iV M on(^=15°) Aft K '' 

and by assuming that the cross section at i?^ = 15° is given by pure Rutherford scattering. An absolute determi- 
nation of the scattering cross section from the accumulated charge and the target thickness agrees within the quoted 
uncertainties with the relative determination using the monitor detectors. 

In the measured angular range the cross section covers more than 4 orders of magnitude. The result of the 
experiment, normalized to the Rutherford cross section, is shown in Fig. [|. The error bars in Fig. || contain statistical 
uncertainties (< 1% for almost any angle) and systematic uncertainties coming from the accuracy of the angular 
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adjustment of the detectors and from contributions of other samarium isotopes in the 144 Sm target (chemical impurities 
of the targets are negligible) . 

At forward angles the accuracy of the angular adjustment (< ±0.1°) leads to an uncertainty of about 1% in the 
determination of the cross section, at backward angles this uncertainty in the cross section practically disappears. In 
contrast, at forward angles the elastic scattering cross section of all samarium isotopes is very close to the Rutherford 
cross section, however at backward angles the cross section measured with the natural samarium target is about 20% 
to 30% smaller than the 144 Sm scattering cross section. Because of the high enrichment (96.52%) of the 144 Sm target 
the resulting uncertainty remains smaller than 1% even at backward angles. Therefore, we renounced of a correction 
of the 144 Sm scattering cross section. 

The resulting high accuracy of this experiment is better than 2% (including both statistical and systematic uncer- 
tainties) even for the very backward angles where the cross section is more than 4 orders of magnitude smaller than 
at the forward angles measured in this experiment. 



III. OPTICAL MODEL ANALYSIS 



The theoretical analysis of the scattering data was performed in the framework of the optical model (OM). The 
complex optical potential is given by 

U(r)=V c (r) + V(r) + iW(r) (3.1) 

where Vc(r) is the Coulomb potential, and V(r) resp. W(r) are the real and the imaginary part of the nuclear 
potential, respectively. 

The real part of the optical potential was calculated by a double-folding procedure: 

Vf(r) = J J p P (r P ) p T (r T )v cS (E,p = p P + p T ,s = \r + rp - rT\)d 3 r P d 3 r T (3.2) 

where pp, p^ are the densities of projectile and target, respectively, and v c g is the effective nucleon-nucleon interaction 
taken in the well-established DDM3Y p aram etrization pi| , p2[ . Details about the folding procedure can be found in 
Refs. 0,0, the folding integral in Eq. ^2 was calculated using the code DFOLD |^3j. The strength of the folding 



potential is adjusted by the usual strength parameter A with A « 1.2 — 1.3. 

The densities of the a particle and the 144 Sm nucleus were derived from the experimentally known charge density 
distributions p4|], assuming identical proton and neutron distributions. For N ~ Z nuclei up to 90 Zr (Z = 40, 
N = 50) this assumption works well, however, in the case of 208 Pb [Z — 82, N = 126) a theoretically derived neutron 
distribution and the experimental proton distribution had to be used to obtain a good description of the elastic 
scattering angular distribution JLl| . To take the possibility into account that the proton and neutron distributions 
are not identical in the nucleus Sm (Z — 62, N = 82) a scaling parameter w for the width of the potential was 
introduced, which is very close to unity. The resulting real part of the optical potential is given by: 

V(r) = A • V{(r/w) (3.3) 

For a comparison of different potentials we use the integral parameters volume integral per interacting nucleon pair 
Jr and the root-mean-square (rms) radius r rmS) R, which are given by: 

Jr = -At~ I V(r)d 3 r (3.4) 

'rms,R — 



A P A T 

-,1/2 

(3.5) 



/ V(r) r 2 d 3 r 



JV{r) d?r 

for the real part of the potential V(r) and corresponding equations hold for W(r). The values for the folding potential 
V { at E uh = 20 MeV (with A = w = 1) are J R = 260.41 MeVfm 3 and r rmS)R = 5.573 fin. 

The Coulomb potential is taken in the usual form of a homogeneously charged sphere where the Coulomb radius 
Rc was chosen identically with the rms radius of the folding potential Vf. Rq = Frms,R = 5.573 fm. 

For the imaginary part of the potential different parametrizations were chosen: the usual Woods-Saxon (WS) 
potential 

W ws (r) = W -{l + exp{r-R)/a)- p (3.6) 
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1 /3 

where R is usually given by R = Rq • Ap , and a series of Fourier-Bessel (FB) functions 

n 

W FB (r) =J2 a k sin (kirr / R FB ) / {knr / R FB ) (3.7) 



fc=i 



with the cutoff radius -Rfb ■ 

A fitting procedure was used to minimize the deviation x 2 between the experimental and the calculated cross 
section: 



i=l 



The calculations were performed using the code AO p5] . 

The parameters of the imaginary part, the potential strength parameter A and the width parameter w of the 
real part, and the absolute value of the angular distribution were adjusted to the experimental data by the fitting 
procedure. The ratio r of the calculated to the experimental absolute value of the angular distribution is very close 
to 1 in all fits: r — 1.006 ±0.002. Therefore, we renormalized the measured angular distribution by this ratio r. This 
renormalization of 0.6% lies well within our experimental uncertainties. 

Four parametrizations of the imaginary part of the potential were used to determine the influence of the type of 
imaginary parametrization on the resulting volume integrals. They are labelled by 1, 2, 3 and 4 in Tabs. | and [y. 
Labels 1 and 2 correspond to WS potentials with p=l and p=2, respectively. For the parametrizations 3 and 4 FB- 
functions were used with 5 and 8 FB coefficients, respectively. The results are listed in Tabs. | and ||, the calculated 
angular distributions are shown in Figs. |§| and ||. In addition, also the angular distribution derived from a "standard" 
WS potential is shown, which was used in previous calculations of the 144 Sm(a,7) 148 Gd capture cross section HQ. 
The fits 1-4 look very similar; the x 2 /F varies from 1.74 to 1.82. However, the fits 3 and 4 using FB functions in the 
imaginary part show a slightly oscillating imaginary potential; fit 4 even shows a small unphysical region where the 
imaginary part becomes positive. 

Of course, the x 2 /F °f fit 5 is inferior compared to fits 1-4 because the WS parameters taken from the study of 
Ref. H were not adjusted to the experimental angular distribution. 



IV. DISCUSSION 



A. Ambiguities of the optical potential 

We extracted a definite optical potential from the experimental elastic scattering data on 144 Sm(a,a) 144 Sm at 
Ei&h = 20 McV. First of all, very accurately measured scattering data are necessary for this determination, and 
furthermore a definite solution has to be selected from several potentials which describe the data almost identically. 
These problems result from discrete ambiguities (the so-called "family problem") and from continuous ambiguities. 

The problem of continuous ambiguities is reduced to a great extent by the use of folding potentials because the 
shape of the folding potential is better fixed compared to standard potentials of WS type. The width parameter of 



the folding potential, which was introduced in Sect. Ill should remain very close to unity; otherwise the parameters 
extracted from this calculation arc not very reliable. 

The "family problem" is illustrated in Fig. |[ Because of the reasons mentioned above we used the simplest 
parametrization of the imaginary potential, the standard WS type with p = 1 as employed in fit 1. If one now varies 
continuously the depth of the real part of the optical potential (i. e. the strength parameter A) and adjusts as well the 
width parameter w and the parameters of the imaginary part to the experimental data, then one obtains a continuous 
variation of w and Jr, but oscillations in x 2 '/F correlated with oscillations in Jj. Each (local) minimum in X 2 IF ~ 
shown as data points in Fig. ^ - corresponds to one family of the optical potential. Obviously, the deepest minima 
in x 2 IF can be obtained for the families 4, 5, and 6 (see Fig. ||). This restriction on the number of the family is 
confirmed by the behaviour of the width parameter w, which should be very close to 1. Fit. 1 in Sec. Ill (sec Figs. |^, 
H) can be found here as family 4. 

A more detailed analysis of the resulting real potentials shows that all potentials have the same depth at the 
radius r — 10.61 fm: V(r — 10.61 fm) = —0.54 MeV. This result is shown in Fig. [| (An exception was found for 
family 1, but for the region of this family x 2 IF does not show a certain minimum, see Fig. [5].) A similar result was 
already found by Badawy et al. [p6[ by analyzing excitation functions of a scattering measured close to d = 180° 
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at energies 10 MeV < £"i a b < 20 MeV: These authors stated that "at energies near the Coulomb barrier the a— 
particle scattering data are fitted by any Woods-Saxon potential whose depth at r = Rq is 0.2 MeV". For 144 Sm 
they derived Rq = 11.04 ± 0.02 fm. However, in contradiction to that statement we point out that one has to 
choose one of the discrete values of the real volume integrals Jr determined by the minima in x 2 1 ' F (families 1-11). 
Choosing e.g. a strength parameter A = 1.35 (which is exactly between families 4 and 5) and adjusting w so that 
V(r — 10.61 fm) = —0.54 MeV (w — 1.018) results in a significantly worse description of the experimental data 
(x 2 1 ' F — 1-92 compared to x 2 /F = 1-82 for families 4 and 5). Of course, this discrimination is only possible if very 
accurately measured scattering data are available! 

The "family problem" usually can be solved at higher energies. The 144 Sm(a,a) 144 Sm scattering data at E\ a ^ = 
120 MeV [^j can be described well by a calculation using a folding potential with Jr = 286.8 MeVfm 3 ]Tl|| , and 
a similar volume integral was found in Ref. [^7| using WS potentials. Together with the systematic behaviour of 
folding potentials for nuclei with A > 90 |ll| which can be described by the interplay of the energy dependence of the 
NN interaction v e f[ and the effect of the so-called "threshold anomaly" fl2],Jl6|] one expects volume integrals of about 



Jfj ps 330 — 340 MeVfm 3 at the low energies analyzed in this work (see Fig. RTupper part and Sect. IVB). /From this 
point of view we can decide that family 4 corresponds to the volume integral Ji?(£q a b = 120 MeV) = 286.8 MeVfm 3 . 

There is a further confirmation for family 4. The ground state wave function of 148 Gd gs . = 144 Sm ® a can be 
calculated fll|| . The number of nodes N and the angular momentum L of the a particle outside the 144 Sm core are 
related to the oscillator quantum number Q by 

4 4 

Q = 2N + L=J2(2ni + k) = J2 qi - ( 4<1 ) 

i=l i=l 

Using Q = 18 (corresponding to qi = 5 oscillator quanta for each neutron in the 2/7/2 shell and qi — 4 oscillator 
quanta for each proton in the 2d 5 / 2 shell, see e.g. |p8[ ) one has to adjust the folding potential strength to reproduce 
the binding energy of the 148 Gd ground state. One obtains A = 1.159 and J R = 311.2 MeVfm 3 at E = +3.2 MeV (the 
nucleus 148 Gd decays by a emission) which fits into the known systematic behaviour of a-nucleus volume integrals 

0- 

For all these reasons we chose the calculations corresponding to family 4 with Jr w 349 MeVfm 3 and Jj « 
52.5 MeVfm 3 (see Sect. HI) obtained with the standard WS parametrization of the imaginary part (fit 1 in Tabs. || 
and||). 



B. Extrapolation to E c . m . = 9.5 MeV 

For the calculation of the 144 Sm(a,7) 148 Gd reaction cross section at the astrophysically relevant energy i? c . m . = 
9.5 MeV one has to determine the optical potential at that energy. The following methods were applied to extract 
the real and imaginary part of the potential. 

In a first step the folding potential in the real part was calcul ated at the energy E c , m , — 9.5 MeV (this is necessary 



because of the energy dependence of the interaction v e s, see Eq. 3.2). Second, the width parameter w was taken from 
fit 1: w = 1.022. Third, because of the rise of the volume integrals Jr at low energies which is about AJr/AE « 
1 — 2 McVfm 3 /MeV jllj we adjusted the parameter A to obtain a volume integral of JR(E c _ m , = 9.5 MeV) = 
334 ± 6 MeVfm 3 : A = 1.1965 ± 0.0216. The uncertainties of Jr and A were estimated from the uncertainties of 
AJr/AE and J R (Eub = 20 MeV). 

The volume integral of the imaginary part can be parametrized according to Brown and Rho (BR) [^9| : 

f for E c m < E Q 

Jl(E c . m .) = \ (E c . m .-Eof , (4.2) 

with the excitation energy Eq of the first excited state, and the saturation parameter Jq and the rise parameter A, 
which are adjusted to the experimentally derived values. 

/From the 144 Sm scattering data at £i a b = 20 MeV and E iah = 120 MeV @ one obtains Jq = 79.98 MeVfm 3 and 
A = 12.84 MeV. The excitation energy of the first excited 2+ state in 144 Sm is Eq = 1.660 MeV. This leads to a 
volume integral at E c . m . = 9.5 MeV of J/(S c .m. = 9.5 MeV) = 21.7 MeVfm 3 . Because of the weak mass dependence 
of the imaginary part for heavy nuclei with a magic neutron or proton number (see Fig. [t|, lower part) and because 
the BR parametrization is not very well defined by 2 data points (and 2 parameters to be adjusted!) we also used the 
well-defined BR parametrization for 90 Zr taken from Ref. Q: Jq = 84.3 MeVfm 3 and A = 11.8 MeV. This leads 
to Ji(E c . m . = 9.5 MeV) = 25.5 ± 0.3 MeVfm 3 using either Eq = 1.66 MeV adjusted to 144 Sm or Eq = 1.78 MeV 
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adjusted to 90 Zr. Combining the values derived from the BR parametrizations of 144 Sm and 90 Zr, we adopt a volume 
integral of Ji(E c , m , — 9.5 MeV) = 22.5^;° MeVfm 3 . 

The shape of the imaginary potential is not as certain as its volume integral. Several parametrizations lead to 
almost identical fits at Ey^ — 20 MeV. For our extrapolation we used the geometry of the potential derived in fit 1 
(see Tabs. | and |l|), because the shape of the imaginary potential should not change dramatically from E\ a h = 20 MeV 
to E c _ m , = 9.5 MeV, and we adjusted the depth Wq = 4-55^030 MeV to obtain a potential with the correct imaginary 
volume integral. 



C. 144 Sm(a,7) 148 Gd and the 146 Sm/ 144 Sm production ratio 



Reaction rates for the reaction 144 Sm(a,7) 148 Gd at Tg = 2.5 and Tg = 3.0 are listed in Tab. III. To determine the 
influence of the optical potential on the 144 Sm(a,7) 148 Gd cross section the calculations of Refs. ||(?J were repeated 
only changing the optical potentials. Compared to a previous calculation using a folding potential with parameters 
derived from a global systematics HQ the reaction rate is reduced by a factor of about 1.5. The optical potential 
is well-defined by the scattering data. With the optical potential determined in this work the calculated and the 
preliminary experimental 144 Sm(a,7) 148 Gd capture cross section agree reasonably well. 

The reaction rate does not depend strongly on the chosen family because the scattering data are reproduced quite 
well from calculations using potentials of families 3, 4, and 5. The capture cross section decreases (increases) by 
about 11% (10%) using family 3 (5) instead of family 4. A similar increase of abo ut 8% compared to family 4 is 



obtained when one uses a potential between families 4 and 5 as mentioned in Sect. IV A because of the somewhat 
larger imaginary part of the potential. 

The influence of the reaction rate on the 146 Sm/ 144 Sm production ratio was analyzed in Ref. 0. Our reduced 
reaction rate lies between cases C and D of Tab. I in Ref. j| corresponding to a production ratio of about 0.3 which 
is well between the experimentally derived limits of 0.1 to 0.7 ||. 



V. CONCLUSION 



The elastic scattering cross section 144 Sm(a,a) 144 Sm was measured at E^ — 20 MeV with very high accuracy. 
A definite optical potential could be derived from the experimental data and by taking into account the systematics 
of a-nucleus optical potentials of Ref. 111]]. Again using this systematics, an extrapolation of the optical potential 
from = 20 MeV to the astrophysically relevant energy -E c . m . = 9.5 MeV was possible with very limited un- 
certainties. The 144 Sm(a,7) 148 Gd cross section is reduced by a factor of about 1.5 compared to a previous folding 
potential calculation, it still lies between the two ESW calculations differing by a factor of 10. The uncertainty of the 
144 Sm(o!,7) 148 Gd cross section coming from continuous and discrete ambiguities of the optical potential is reduced 
by a great amount. The use of systematic folding potentials is highly recommended for the analysis of low-energy 
a scattering and a capture reactions because of the reduced number of free parameters compared to previous calcu- 
lations using WS potentials, resp. because of the reduced uncertainties in the radius parameter compared to ESW 
calculations. 
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TABLE I. Potential parameters of the imaginary part of the optical potential derived from the angular distribution of 
4 Sm(a,a) 144 Sm at £i ab = 20 MeV. 



fit 


W (MeV) 


#o (fm) 


a (fm) 


P 










1 


10.64 


1.6758 


0.1680 


1 










2 


10.70 


1.7132 


0.2265 


2 










fit 


Rfb (fm) 


ai 


0.2 


a3 








as 


3 


15.0 


-11.65 


-16.66 


6.09 


21.31 


10.79 






4 


15.0 


-9.50 


-6.50 


19.42 


14.84 


-23.51 


-35.83 -14.88 


-1.06 


5 a 




V = 


185 MeV, W 


= 25 MeV, Ro, R = Ro,i 


= 1.4 fm, a_R 


= a 7 = 0.52 fm 





Ref. § 



TABLE II. Integral potential parameters J and r rms of the real and imaginary part of the optical potential derived from 
the angular distribution of 144 Sm(a,a) 144 Sm at E uh = 20 MeV. 



fit 


A 




Jr 


^rms,R 


Ji 




X 2 /F 


No. 






(MeVfm 3 ) 


(fm) 


(MeVfm 3 ) 


(fm) 




1 


1.2568 


1.0220 


349.34 


5.6961 


52.57 


6.8311 


1.823 


2 


1.2580 


1.0216 


349.29 


5.6940 


52.41 


6.8142 


1.823 


3 


1.3425 


0.9976 


347.04 


5.5600 


57.49 


6.0756 


1.807 


4 


1.2771 


1.0067 


339.31 


5.6110 


51.30 


6.8868 


1.738 


5 a 


Woods- 


-Saxon 


557.59 


6.0026 


75.35 


6.0026 


41.0 



Ref. | 



TABLE III. Reaction rates of 144 Sm(a,7) 148 Gd at temperatures of T 9 = 2.5 and T 9 = 3.0. 



potential 


potential 


(«,y) 


N A - < <7V > 


(cm s mole ) 




from Ref. 


from Ref. 


T 9 = 2.5 


T 9 = 3.0 

r, ro . . 1 r.-ia 



ESW, R = 8.75 fm 
WS 

Folding, A = 1.159 
Folding 



m 

this work 



this work 



3.75 x 10" 15 

1.95 x 10~ 15 

1.27 x 10" 15 

7.91 x 10" 16 



2.35 x 10~ 12 

1.22 x 10~ 12 

7.56 x 10 -13 

5.63 x 10~ 13 







2.5 



2.0 



2 1-5 

CD 
• i— i 



% 1.0 



0.5 



0.0 



_L 



-o 
o 



d 

DC 



o 

O 
<U 
O. 

X 



41.0 41.5 42.0 42.5 

^iab (deg) 



43.0 



43.5 



FIG. 1. Relative yield of 1 C recoil nuclei in coincidence with elastically scattered a particles. The shaded area shows the 
angle and the uncertainty which is expected from the calibration using the steep kinematics of 1 H(a,a) 1 H. The dotted line is 
a Gaussian fit to the experimental data points to guide the eye. 
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FIG. 2. Typical spectra of 144 Sm(a,a) 144 Sm at $i a b = 25° (upper diagram) and at $i a b = 130° (lower diagram). 
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FIG. 3. Elastic scattering cross section of 144 Sm(a,a) 144 Sm normalized to the Rutherford cross section. The line is the 
result of an OM calculation using folding potentials and corresponds to fit 1 of Tabs. 1 and 2 (see Sect. III). The inserts show 
magnifications of the forward and backward angular range. 
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FIG. 4. Calculated elastic scattering cross section of 144 Sm(a,a) 144 Sm normalized to the Rutherford cross section. The 
lines are the result of different OM calculations using folding potentials (fits 1-4, which have almost the same / F) and using 
a standard WS potential (fit 5) from Ref. M. The parameters are shown in Tabs. 1 and 2. 
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FIG. 5. Parameters of different families (numbered 1-11) of optical potentials derived from the 144 Sm(a,a) 144 Sm scattering 
data: strength parameter A, width parameter w, volume integrals of the real and imaginary part Jr resp. Ji, and the deviation 
per degree of freedom x 2 / F (from top to bottom). The data points are results of fits of real and imaginary part using 
different starting values, the lines are the result of an interpolation in A and w to adjust the real part of the potential to 
V(r = 10.61 fm) = —0.54 McV, and the imaginary part of the potential was again adjusted to the experimental data. The 
deepest minima in x 2 /F are obtained for the families 4, 5, and 6. 
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FIG. 6. Real part of different families of optical potentials (labelled 1-11) derived from the 144 Sm(a,a) 144 Sm scattering data: 
in the insert one can see that all potentials have the same depth at the radius r — 10.61 fm: V(r = 10.61 fm) = —0.54 MeV. 
For family 1 (dashed line) one finds an exception but for the region of this family x 2 / F does not show a certain minimum (see 
Fig.|). " 
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FIG. 7. Volume integrals of the real (upper) and imaginary part (lower diagram) of the optical potential derived from 
144 Sm(a,a) 144 Sm scattering. For comparison the volume integrals derived from 90 Zr(a,a) 90 Zr scattering jll| were added. The 
lines in the lower diagram show the results of BR parametrizations of the imaginary part (see text). 
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